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f~| ■ Abstract. We propose a feasible scheme for generation of strongly non-Gaussian 

states using the cross-Kerr nonlinearity. The resultant states are highly non-classical 
^ ' states of electromagnetic field and exhibit negativity of their Wigner function, sub- 

. Poissonian photon statistics, and amplitude squeezing. Furthermore, the Wigner 

^ | function has a distinctly pronounced "banana" or "crescent" shape specific for the 

Kerr-type interactions, which so far was not demonstrated experimentally. We show 
that creating and detecting such states should be possible with the present technology 
using electromagnetically induced transparency in a four-level atomic system in re- 
configuration. 



PACS numbers: 42.65.-k, 42.50.-p, 03.67.-a 
1. Introduction 



Non-Gaussian and in general non-classical states of light enjoy an increasing interest 
in the quantum information community. They are indispensable building blocks 
for entanglement generation, distillation and broadcasting, for quantum information 
relays and networks in the realm of quantum information processing based on infinite- 
dimensional optical systems and atomic ensembles [H [2]. As has been shown, it is 
impossible to perform optical quantum computation with linear optics only and a non- 
linear interaction is needed to perform a universal quantum computation pJJ. Hence 
it is unavoidable to investigate the nonlinearities of atomic systems or alternative 
solutions for giant nonlinearities and efficient nonlinear coupling, such as measurement- 
induced nonlinearities or other non-Gaussian operations. In addition, the experimental 
difficulty associated with the "on-line" nonlinear operations can be circumvent by using 
the off-line resources, the highly non-Gaussian states prepared separately from the 
actual computation and fed into the computation circuits when necessary to replace the 
nonlinear gates [U [5] . Continuous- variable quantum repeaters, quantum relays and a 
number of other computation and communication tasks rely on non-Gaussian operations 
or non-Gaussian states, states described by a non-Gaussian Wigner function. 
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One of the examples of a non-linear coupling is the cross-Kerr effect which involves 
controlling the refractive index experienced by one mode of the electromagnetic field by 
the intensity of another. In this paper, we propose an experiment that employs the cross- 
Kerr effect to create highly non-classical non-Gaussian states of light via interaction 
of two coherent beams in an atomic medium exhibiting electromagnetically-induced 
transparency [6], subsequent measurement on one beam and feed- forward on the other. 
Such experiment on the one hand responses to the need for an efficient source of non- 
Gaussian states, on the other hand can be seen as a test-bench for the strong non-linear 
coupling using an EIT based system. In our case, the first evidence of the large non-linear 
phase shift in combination with enough coherence to generate and preserve quantum 
features can be obtained using merely a direct photodetection to verify the photon- 
number squeezing in one of the output modes. This would pave the way for further 
applications of such non-linear systems: Cross-Kerr nonlinear interaction provides a 
basis for several proposals of quantum information protocols or their elements, such as 
non-demolition photon number detection |7|, C-NOT gate [8J, or continuous- variable 
entanglement concentration [9| [TP]. 

The first basic ingredient of our scheme is the nonlinear coupling based on the 
third-order nonlinearity, an optical cross-Kerr effect [11]. Such Kerr-type interaction 
of the two initially independent coherent beams a and b entangles them producing a 
continuous-variable state of modes a and b with quantum correlations between photon 
number in one mode and phase of the other. A subsequent measurement done on the 
mode b generates a conditional squeezed state of the mode a. This is a photon-number 
squeezed state described for the first time in [12]. The mechanism behind it is a re- 
shaping of the quantum uncertainty due to the Kerr-effect and the influence of a local 
measurement on the system of two spatially-separated entangled modes. A particular 
feature of the states squeezed using Kerr nonlinearity is their non-Gaussian character in 
contrast to the squeezed states produced, e.g., in the parametric process based on the 
nonlinearity of the second order [2]. The corresponding Wigner function has a specific 
crescent shape and exhibits negativity in some regions of the phase space in the form 
of the decaying fringes, resembling a section of the Wigner function of the Fock states. 
Indeed, there is a strong connection between the states emerging in this scheme and the 
Fock states: the generation of the non-classical state in our scheme can be understood 
as a superposition effect between different photon-number or Fock states. Hence the 
resultant Wigner function can be viewed a result of interference between the Fock-state 
Wigner functions with different photon numbers n. We will provide a more detailed 
discussion of the mechanism producing this particular state in Sec. 3 and 4 of the 
paper. 

Another inherent feature of our protocol is the controlled displacement. We are 
dealing with continuous variable quantum systems. It means, that although in each 
single run of the proposed experiment (corresponding to a single measurement result 
on mode b) a squeezed state is created, the detection of squeezing will not be possible 
because it requires a statistical processing of an ensemble of measurement results. In 
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a number of subsequent runs, the overall squeezing vanishes, as we then will obtain a 
mixture of quantum uncertainties with different mean photon numbers and the overall 
state will become fuzzy. The way out is provided by the quantum feed-forward: the state 
in the mode a is displaced in each run with a displacement amplitude determined by a 
corresponding measurement outcome in mode h. This procedure merges all differently 
displaced crescent-shaped uncertainties into a single crescent which exhibits amplitude 
(or photon-number) squeezing. 

The non-classical non-Gaussian states of this type have not yet been demonstrated 
experimentally although the Kerr effect, e.g. in optical nonlinear fibres, was exploited 
extensively for squeezing and entanglement generation (see e.g. [EH IE])- The 
challenging point is to produce the strong enough nonlinear coupling to generate sensible 
non-Gaussian features. All the experiments performed so far were working in the regime 
of weak nonlinearity restricting themselves to the first stages of quantum state evolution, 
which can be well approximated by the Gaussian Wigner function of the same type as 
the one describing squeezing in parametric second-order nonlinear processes. The most 
promising candidate for observing the large cross-Kerr effect is the four-level atomic 
medium exhibiting electromagnetically induced transparency (EIT) [6]. Under certain 
conditions, the cross-phase modulation in such media can reach values by many orders 
of magnitude larger than it is possible in optical fibres and the phase shift between 
two photons can reach values of 10~ 3 for hot atoms in a vapour cell or even 10 _1 for 
cold atoms in a magneto-optical trap [15] . The experimental feasibility of achieving 
a large cross-Kerr interaction for continuous- wave fields has been demonstrated using 
cold atoms in magneto-optical trap [IB] . Recently, a suitable method was suggested for 
group- velocity matched pulses and room temperature atomic gas cell [15]. The form of 
the cross-Kerr interaction in a such four-level atomic system in the N-configuration can 
be derived in a simple way [17]. Essentially, the strong cross-Kerr interaction is caused 
by the ac-Stark shift produced by the third perturbing field to the dark-state of the 
lambda subsystem. 

The paper is organised as follows. In Sec. [2] we present the experimental scheme. 
Sec. |3] discusses the non-Gaussian properties of the state obtained in a single run of the 
experiment. In Sec. H]we show how a controlled displacement on this output state can 
reproduce the same highly non-Gaussian state in multiple runs generating a constant 
output state with high accuracy. In Sec. [5] we discuss the scaling of the non-Gaussian 
effects with experimental parameters and conclude in Sec. [6j 

2. Experimental scheme 

Consider the following experimental scheme sketched in Figure [TJ two coherent states 
with amplitudes a and (3 occupying modes a and b, respectively, interact in a cross- 
Kerr medium and mode b is then subject to a measurement of the x quadrature via 
homodyne detection. Based on the measurement outcome x, a displacement operation 
is performed on mode a. In this section we obtain the form of the output state p and 
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Figure 1. Proposed experimental setup for generating crescent states with negative 
Wigner function: two coherent states interact in a cross-Kerr medium, one of them is 
subject to measurement of the field quadrature, and displacement is performed on the 
other one depending on the measurement outcome. 
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its properties will be discussed in the following sections. 
The Hamiltonian of the cross-Kerr interaction is 

H = -ih a h bl (1) 

where 7' expresses the strength of the interaction and h a = d^a, h b = b^b are the 
photon number operators of the two modes. We introduce the phase shift between two 
photons as 7 = j't, where t is the interaction time, which gives the evolution operator 
U = exp(ijn a n b ) . If the two modes are originally in coherent states \a) a and \/3) b , then 
after the cross-Kerr interaction their state will be 

2 °° a n 

\^)ab = exp(jh a h b )\a) a \l3) b = e~ |a| /2 V —j= exp(~fh a h b )\n) a \(3) b 

n=0 

00 n 

= e-N 2 / 2 V4=|n) a |/3e^) 6 (2) 
n V n\ 

This is an entangled state where the photon number n a in mode a is correlated with the 
amplitude of coherent state in mode b by introducing to it a phase shift dependent on 
n a (and vice versa). 

In the next step the quadrature x = (b+&) / y/2 of mode b is measured via homodyne 
detection with the outcome x, which results in the following (unnormalized) state of 
mode a: 

mx)) . = tm + - ^ t °' bpH ' " ^ + ^ " lM ' J I")- (3) 



n=0 v^! 



where we have denoted (3 n = Re (j3e nn ), (3' n = Im (j3e lin ) and used the expression for a 
coherent state in x- representation. The square of the norm of \ip{x)) a is the probability 
density P(x) that the particular measurement outcome x occurs. 

As the last step, a displacement in the phase space is applied to mode a. The 
resulting state of mode a is then 

= D[d(a;)]|^(aO) = e d{x) ^- d ' {x) ~ a \^{x)) , (4) 

where d(x) is the displacement parameter that depends on the measurement outcome 
x. When the experiment is performed repeatedly, then the output state is averaged to 

p= I m x ))(<S>(x)\dx, (5) 
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which is properly normalized. 

The state in general depends on the measurement outcome x. However, we 

will show that in a certain regime this dependence can be quite weak. This way, by 
running the experiment many times (and obtaining different values of x), one can still 
get almost a pure state p at the output mode a that has highly non-classical properties. 

3. Properties of the state \4>(x)) 

Before describing the averaged output state p, in this section we will focus on the state 
\i/j(x)) of mode a after the x measurement on mode b has been performed. Consider 
the specific situation of (3 — i|/3|e~ 17 ' a ' 2 , i.e., when the phase of of the mode b is set to 
7r/2 — 7|a| 2 . We also put a constraint on \a\, such that \a\ < 10 and the product |a/3|7 
is of order of unity. The reasons for these assumptions will be explained later in this 
section. Taking into account that 7 is several orders of magnitude smaller than unity, 
it also holds that j\a\ <C 1. 
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Figure 2. Photon number probability distribution of the coherent state \a) (a, 
black curve) compared with the normalized distributions of the states \%p(x)). The 
parameters are \a\ = 6, j\/3\ = 0.36 and x = —4 (b, green curve), x = (c, blue curve) 
and x = 4 (d, red curve). The photon number distribution in the states \tp(x)} is about 
4x squeezed compared to the Poissonian statistics of the coherent state |a). 

In this particular case, for evaluating the real and imaginary parts of /3e 17n to 
calculate (3 n and f3' n in ([3]), we use the expansion of the exponential function as follows: 

/3e i7n = i|/3|e i7(n - |a|2) w i|/3| - 7 |/3|(n - \a\ 2 ). (6) 

We terminated the series after the second term since 7(71 — \ct\ 2 ) 1 for all n for which 
the probability of having n photons in the state |a) is non-negligible. Indeed, the photon 
number distribution of the state \a) is Poissonian with both mean and variance equal to 
|a| 2 , hence (n — \a\ 2 ) is of order of |a| and therefore j(n — \a\ 2 ) 1 holds in agreement 
with our assumption j\a\ <^ 1. 

From (jHJ) we get f3 n = — |/3|7(n — |a| 2 ) and f3' n = \/3\. Substituting this into (j3J) we 
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Figure 3. Illustration in the phase space of how the photon number squeezing of the 
state \ip(x)} emerges: the red circles represent the uncertainty regions of coherent states 
|e 17 "/3)b in ([2|) entangled with the photon numbers n in mode a (n fa \a\ 2 denotes the 
most likely photon number). The circles are distributed along the green arch centred 
at the origin of the phase space, having the radius \j3\ and stretching to the angle of 
approximately -7 |ct| to both sides of the imaginary axis. If then the quadrature x is 
measured on mode b with a particular outcome x (blue line represents the quadrature 
eigenstate \x}), then only those states are chosen from the superposition ([3|) that have 
a non-negligible overlap with |a;). If the length of the arch, 27|a/?|, is larger than 
the spatial extension of the uncertainty region of about unity, then only a relatively 
small number of circles have a significant overlap with \x) and so the photon number 
variance in the state is reduced significantly compared to the state \a). The 

effect is the strongest when the arch is, very loosely speaking, perpendicular to the line 
representing |x); this gives the phase condition on 0. 
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\n), (7) 



where <\> = y2x\/3\ + 7|a/3| 2 is an irrelevant global phase factor. The important term 
here is the exponent, which provides the insight into the main physical effects due the 
cross-Kerr interaction. The first term in the exponent shows that the phase of the state 
\ijj(x)) is increased by 7|/3| 2 with respect to the phase of the original state \a). The 
second, Gaussian term causes the photon number distribution in \ip(x)) to be altered 
with respect to the original Poissonian distribution of the state \a): If 2|a/5|7 > 1, then 
this term becomes more important than the factor a n /\fn\. Moreover, as the product 
2|a/3|7 grows above unity, the photon number distribution quickly gets dominated by 
the Gaussian term. In this situation the mean photon number is \a\ 2 — x/{\/2\(5\^) 
compared to | or | 2 for the coherent state |a) and the variance of the photon number 
distribution is l/{2^\f3\) 2 compared to \a\ 2 for \a). This means that the state \ip(x)) 
is squeezed in the photon number by the factor of approximately 2|a/3|7, as illustrated 
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Figure 4. Two different views of the Wigner function W(u) of the output state \if>(x)) 
for \a\ = 6, -f\f3\ = 0.36 and x — 0. The product \a/3j\ is equal to 2.16 and the phase 
of a is set to — j\(3\ 2 to make the phase of the output state zero. The crescent shape 
and negative values of the Wigner function are clearly visible. 

in Figure [2] for |a/3|7 ~ 2.16 and different values of the measurement result x. Thus 
if 2|a/3|7 > 1, then the output state will exhibit a strongly sub-Poissonian photon 
statistics. Figure [3] gives a visualisation of the above argumentation. Equation ([7]) 
and Figure [3] also provide an argument for the particular choice of the phase of the 
mode b, (3 = i\p\e-^ a \ 2 . This way we have demonstrated the photon number variance 
reduction - photon number squeezing in the output state \ip{x)) obtained in a single run 
of the experiment. A consequence of this reduction are several non-classical phenomena: 
negativity and crescent shape of the Wigner function, as well as amplitude squeezing. 

Figures H] and [5] show the Wigner function corresponding to the output state \ip{x)) 
of ([7j) for different values of \a\ and j\/3\. The above mentioned features of the Wigner 
function are clearly visible especially in FigureHl Comparison of FiguresH]and[5]provides 
a clear insight into the assumptions made in the beginning of this section. For large 
values of |a| (say over 10) one gets a state that is close to a squeezed Gaussian state 
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Figure 5. Two different views of the Wigner function W(u) of the output state \if>(x)) 
for \a\ — 30, -f\f3\ = 0.066 and x = such that \a(3j\ = 1.98. The state is strongly 
squeezed and is close to a Gaussian state, hence the negativeness of the Wigner function 
is weak. 



and the features we are looking for, namely the negativity and crescent-shape of the 
Wigner function, are suppressed (Figure [5]). We also assumed that the product |a/?|7 
is of order of unity because otherwise the photon number squeezing is small, the state 
\i/j(x)) is close to a coherent state and hence does not have strong non-Gaussian and non- 
classical properties. |a/3|7 > 1 requires a large cross-Kerr nonlinearity and/or a large 
amplitude (3. The former we anticipate to be feasible with the four-level EIT based 
system mentioned in the Introduction. The Wigner function in Figure H] is plotted for 
the parameters satisfying all the discussed restrictions, with |a/?|7 ~ 2.16 as in Figure[2j 
Such features as the negativity and crescent-shape are highly pronounced, confirming 
the strongly non-Gaussian character of the output state. 
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4. Effect of displacement in the phase space 

As we have seen, the output state \ip(x)) depends on the value x obtained by the 
measurement. Specifically, for different x, there will be different mean photon numbers 
in the output state of the mode a. To prepare an identical crescent state repeatedly, 
one strategy would be to post-select the output at some sufficiently narrow interval of x, 
which would, however, reduce the success rate significantly Another option, however, 
is to correct for the change of the mean photon number caused by the x-measurement. 
It turns out that this correction can be achieved with a good accuracy by performing a 
phase-space displacement operation D[(d(x)] with displacement parameter d depending 
on x. In this way one can obtain a state almost independent of the measurement 

outcome x, which makes p in ([5]), Figure [1] nearly a pure state. 

To find the displacement magnitude \d(x) \ that would compensate for the change of 
the mean photon number, we recall that for coherent states, the mean photon number 
depends quadratically on the amplitude. The state \if>(x)) has the mean photon number 
| o; | 2 — x/(y2\/3\ r y) and therefore its mean amplitude \a'\ is equal to the square root of 
this number. In order to revert to the original amplitude |a|, we need to displace by 
| a | — \a'\. The direction of the displacement in the phase space is given by the phase of 
the state \ip(x)), which is, as seen from ([7]), equal to 7|/3| 2 + arg(o). This way we arrive 
at the displacement parameter 



To see that the state |$(a;)) in (jlj) is really almost independent of the measurement 
outcome x, we investigate the normalized scalar product 

F(x) = WQ)I»W> (9) 

V<t(0)|t(0)>{4.(x)|4.(x)> 

that expresses the overlap of the displaced output state corresponding to an arbitrary 
x and the one corresponding to x — 0. Figure [6] shows the function F(x) and the 
probability density P(x) = (ip(x)\ijj(x)) of the quadrature for a few values of |a| and 
j\/3\. It can be seen that in the regions of x where the probability P(x) is non- negligible, 
the function F(x) has values close to unity, and this is a typical behaviour. This means 
that the state |$(a;)) is close to |$(0)) for all x that are likely to be found in the 
quadrature measurement. Therefore p, which is the mixture of these states [see flSJ)], 
will be almost a pure state. Numerical simulations confirm this and the purity V = Tr p 2 
has values of approximately 0.95 for the states corresponding to parameters a, (3 and 7 
in Figures H] and [51 

Experimentally, the desired controlled displacement can be realized mixing the state 
\i/j(x)) with a coherent state \d/t) on a beam splitter with a very low transmissivity t. If 
\t\ <C 1, then one obtains almost exactly the state D(d)\ip(x)) at one beam splitter 
output port [18]. The amplitude d/t can be varied electronically according to the 
measured value x and (JSj) using, e.g., the scheme depicted in Figured 
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Figure 6. The normalized scalar product F(x) (solid red line) and the probability 
density P{x) (dashed green line) as a function of x for (a) \a\ — 6,j\/3\ — 0.36 and 
(b) \a\ — 9,"f\f3\ — 0.2. Since the function F(x) is much broader than P(x), the 
fidelity F(x) is almost unity for all values x likely to be found in the measurement, 
and therefore the state is almost independent of x. 




Figure 7. A possible way of realizing the displacement operation: a vertically 
polarised coherent beam is directed into a Pockels cell (PC) with the optical axis 
rotated by 45 degrees and controlled by the outcome x of quadrature measurement. 
The output beam then impinges onto a polarising beam splitter (PBS), which yields 
a polarised coherent beam with amplitude depending on x. This beam is then mixed 
with the state \ip(x)) on a highly reflective beam splitter (BS). The displaced state 
\&(x)} is found at one output. 



5. Scaling of the non-classical effects with |a|, \/3\ and 7 

In what follows, we discuss the scaling of the effects described above for the output state 
|$(0)) = \ip(0)), as we have shown that the state is almost independent of x. An 

inspection of formulas (JTj) and © shows that, apart from phase factors, the form of the 
output states depends on |a| and the product "/\f3\. 

Consider first the effect of varying 7|/3| while keeping \a\ constant. For j\{3\ = 0, the 
output state reduces to the original state \a). For very small j\/3\ for which j\a/3\ <C 1, 
the state \ip(0)) is close to a weakly squeezed coherent state with amplitude a exp(ij\(3\ 2 ), 
so it is almost Gaussian. As 7|/3| increases, |^(0)) starts to deviate from a Gaussian state, 
negative regions of the Wigner functions get larger, and the crescent shape of \ip(0)) 
emerges, as well as the photon number squeezing. When 7|/3| gets larger than unity, 
\ip{x)) is almost a Fock state \n) with photon number n depending on the measured 
value of x, and it is no longer true that |$(x)) is almost independent of x. In this 
situation our scheme (without the displacement operation) can be used as a conditional 
source of Fock states. 

The effect of varying a while keeping the product 7|a/5| constant is different. The 
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phase of a influences just the phase of ^(O)). If |a| is small, the crescent shape of the 
Wigner function is strongly pronounced, as well as its negativeness, because the state 
is squeezed in photon number and is not far from a Fock state. Close to the origin of 
the phase space the circles corresponding to a fixed photon number are more curved 
than those farther from the origin, and hence the Wigner function corresponding to 
smaller \oc\ shows a stronger crescent shape than the one corresponding to larger |a|. 
In contrast, for larger |a| the state is closer to a Gaussian state whose Wigner function 
is positive. Hence the negativeness of the Wigner function is more profound in case of 
smaller |a|, and this case is more appealing experimentally. Figures H]- [5] illustrate this 
behaviour. 

6. Conclusion 

In conclusion, we have suggested a feasible scheme to generate highly non-Gaussian 
states exhibiting crescent-shaped Wigner function with negative regions. We envisage 
the application of such states in quantum information processing using infinitely- 
dimensional, continuous-variable quantum systems. The non-Gaussian operations and 
non-Gaussian states attract currently an increasing attention in this context (see e.g. 
|H [19] and references therein) and there is a high demand in the community for their 
experimentally viable implementation. In our scheme, the initial states as well as 
the detection (Figure [I]), feed- forward (Figure [7]) and evaluation steps belong to the 
standard toolbox of quantum optics and are readily available in the laboratory. The first 
verification of quantum features can be performed with merely direct photodetection to 
prove the photon-number squeezing of Fig [2j The negativity of the Wigner function 
and its crescent shape can be visualised using quantum state tomography [20], a 
more involved but established procedure [211 [22]. For the first demonstration of the 
effect using the tomographic measurement of the Wigner function of the state p, the 
displacement operation of Figure can be simplified via replacing it by the electronic 
one. 

The challenging part of the proposed scheme is the strong nonlinear coupling 
implied. Again, it represents an important building block in quantum information 
processing attracting recently a substantial interest from both the theoretical and 
experimental sides. A feasible nonlinear coupling device would make a profound impact 
on the development of quantum communication and computation protocols. So far, 
there are only few first quite involved implementations demonstrating this effect [16J, 
|23j . For the realization of our scheme, we suggest the nonlinear optical cross-Kerr 
effect in an EIT based four-level atomic system, which was closely studied recently 
both theoretically P, US] El El] and experimentally [IB] . Furthermore, due to the 
experimental availability of the other elements and particular simplicity of the first 
verification using direct photodetection, our suggested scheme can serve as a test bench 
for the strong nonlinear coupling preserving quantum effects. 
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